
Calhoun 

iniQiuiic^iul Ar{hiv« of tilt Mil vdl Poii^roduiit School 


Calhoun: The NPS Institutional Archive 
□Space Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


2004-09 

Optimal trajectory reconfiguration and 
retargeting for the X-33 reusable launch vehicle 


Shaffer, Patrick]. 

Monterey California. Naval Postgraduate School 


http://hdl.handle.net/10945/1347 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 

Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


htt p://w ww. n ps. e du/l ib ra ry 


Caflwuo is the Naval Postgraduate School's public access digital repository for 
research mate rials and institutiional publicatkios created by the NPS community. 
Calhoun is named for Professor of Mathematics Guy K. Caftiouo, NPS's first 
appointed — and published — schoteily author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 Univefsity Circle 
Monterey, California USA 93943 







NAVAL 

POSTGRADUATE 

SCHOOL 

MONTEREY, CALIEORNIA 


THESIS 


OPTIMAL TRAJECTORY RECONFIGURATION AND 

RETARGETING FOR THE X-33 REUSABLE 

LAUNCH VEHICLE 


by 


Patrick J. Shaffer 


September 2004 


Thesis Advisor: 

I. Michael. Ross 

Seeond Reader: 

David B. Doman 


Approved for public release; distribution is unlimited 




THIS PAGE INTENTIONALLY LEET BLANK 



REPORT DOCUMENTATION PAGE 


Form Approved OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including 
the time for reviewing instruction, searching existing data sources, gathering and maintaining the data needed, and 
completing and reviewing the collection of information. Send comments regarding this burden estimate or any 
other aspect of this collection of information, including suggestions for reducing this burden, to Washington 
headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 
1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project 
(0704-0188) Washington DC 20503. 


1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

September 2004 Master’s Thesis 


4. TITLE AND SUBTITLE: 5. EUNDING NUMBERS 

Optimal Trajectory Reconfiguration And Retargeting For The X-33 

Reusable Launch Vehicle. _ 

6. AUTHOR(S) Patrick J. Shaffer, CDR, USN 


7. PEREORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Naval Postgraduate School 

Monterey, CA 93943-5000_ 


9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

AFRL/VACA 

Wright Patterson AFB, OH 


8. PERFORMING 
ORGANIZATION REPORT 
NUMBER 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the official 
policy or position of the Department of Defense or the U.S. Government. 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 

Approved for public release; distribution is unlimited 


13. ABSTRACT (maximum 200 words) 

This thesis considers the problem of generating optimal entry trajectories for a reusable launch vehicle following a 
control surface failure. The thesis builds upon the work of Dr. David Doman, Dr. Michael Oppenheimer and Dr. 
Michael Bolender of the Air Vehicles Directorate, Air Force Research Lab Dayton Ohio. The primary focus of 
this work is to demonstrate the feasibility of inner loop reconfiguration and outer loop trajectory retargeting and 
replanning for the X-33 reusable launch vehicle (RLV) following the imposition of a control surface failure. The 
trajectory generation model employs path constraints generated by an AFRL trim deficiency algorithm coupled 
with an inner loop control allocator and aerodynamic database that captures the full 6-DOF vehicle aerodynamic 
effects while utilizing an outer loop 3-DOF model. The resulting optimal trajectory does not violate the trim 
deficiency constraints and provides additional margins for trajectories flown during failure conditions. The 
footprints generated by the thesis show that contemporary footprint analysis for vehicles experiencing control 
surface failures are overly optimistic when compared to those footprints that consider vehicle aerodynamic 
stability and realistic laudable attitudes at the threshold of the landing runway. The results of the thesis also show 
the performance reductions resulting from decoupling the inner and outer loop and that trajectories can be 
generated to the landing runway without using a region of terminal area energy management. 


14. SUBJECT TERMS Reentry, X-33, Trajectory, Optimal, Replanning, Retargeting 


17. SECURITY 
CLASSIFICATION OF 
REPORT 

Unclassified 


NSN 7540-01-280-5500 
Prescribed by ANSI Std. 239-18 


18. SECURITY 
CLASSIFICATION OF THIS 
PAGE 

Unclassified 


19. SECURITY 
CLASSIFICATION OF 
ABSTRACT 

Unclassified 


15. NUMBER OF 
PAGES 

151 


16. PRICE CODE 


20. LIMITATION 
OF ABSTRACT 


Standard Form 298 (Rev. 2-89) 


1 




























THIS PAGE INTENTIONALLY LEET BLANK 


11 



Approved for public release; distribution is unlimited 


OPTIMAL TRAJECTORY RECONFIGURATION AND RETARGETING FOR 
THE X-33 REUSABLE LAUNCH VEHICLE 


Patrick J. Shaffer 
Commander, United States Navy 
B.S., Old Dominion University, 1988 


Submitted in partial fulfillment of the 
requirements for the degrees of 


ASTRONAUTICAL ENGINEER 
AND 

MASTER OF SCIENCE IN ASTRONAUTICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
September 2004 


Author: Patrick J. Shaffer 


Approved by: I. M. Ross 

Thesis Advisor 


David B. Doman 
Second Reader/Co-Advisor 


Anthony J. Healey 

Chairman, Department of Mechanical and Astronautical 
Engineering 



THIS PAGE INTENTIONALLY LEET BLANK 


IV 



ABSTRACT 


This thesis considers the problem of generating optimal entry trajectories for a 
reusable launch vehicle following control surface failures. The thesis builds upon the 
work of Dr. David Doman, Dr. Michael Oppenheimer and Dr. Michael Bolender of the 
Air Vehicles Directorate, Air Force Research Lab Dayton Ohio. The primary focus of 
this work is to demonstrate the feasibility of inner loop reconfiguration and outer loop 
trajectory retargeting and replanning for the X-33 reusable launch vehicle (RLV) 
following the imposition of a control surface failure. The trajectory generation model 
employs path constraints generated by an AFRL trim deficiency algorithm coupled with 
an inner loop control allocator and aerodynamic database that captures the full 6-DOF 
vehicle aerodynamic effects while utilizing an outer loop 3-DOF model. The resulting 
optimal trajectory does not violate the trim deficiency constraints and provides additional 
margins for trajectories flown during failure conditions. The footprints generated by the 
thesis show that contemporary footprint analysis for vehicles experiencing control surface 
failures are overly optimistic when compared to those footprints that consider vehicle 
aerodynamic stability and realistic landable attitudes at the threshold of the landing 
runway. The results of the thesis also show the performance reductions resulting from 
decoupling the inner and outer loop and that trajectories can be generated to the landing 
runway without using a region of terminal area energy management. 
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I. INTRODUCTION 


A. HISTORICAL BACKGROUND 

As early as 1931, Dr. Eugen Sanger eoneeived of an aerovehiele that eould 
aehieve low-earth orbit (LEO) and return to earth like a glider. The Sanger Silverbird 
was advoeated as a hypersonie vehiele that eould deliver an 8300 pound payload into 
LEO and return to earth in an unpowered glide. Alternatively, the Silverbird (or its 
military inearnation, the Amerika Bomber) eould deploy a 1600-pound payload into a 
suborbital trajeetory from a launeh point halfway around the globe and return to its 
origination by skip-gliding or dynamie soaring [ref 25]. 

The problem of extreme reentry heating was severely underestimated by Sanger 
and eonelusions from NACA Ames reentry studies in the 1950’s led to the abandonment 
of skip gliding or dynamie soaring beeause there was no eontemporary teehnologieal 
solution to mitigate the severe heating during reentry. Later manned missions sueh as 
Apollo overeame the reentry heating problem by employing a vehiele that used an 
ablative heat shield to proteet the payload and followed a steep ballistie reentry trajeetory. 
Unfortunately this type of teehnology, while effeetive, did not allow for reusable 
vehieles. 

Interest in reusable spaeeeraft resulted in the development of the spaee shuttle 
whieh makes use of advaneed materials and eontrols teehnologies to fly trajeetories 
whieh keep the shuttle within a very narrow reentry eorridor to limit the heating and 
dynamie pressure effeets of the entry trajeetory thereby ensuring a sueeessful landing 
[ref. 22]. These reentry trajeetories are eomputed offline for a variety of initial 
eonditions and then stored in the shuttle’s eomputer where they ean then be aeeessed by 
the guidanee eomputer. The generation of these trajeetories represents a signifieant effort 
and is supported by a large body of work from eontemporary optimal eontrol researeh. 

Reusable launeh vehieles are being pursued as a low-eost alternative to 
expendable launeh vehieles and the shuttle. The eoneept of a reusable launeh vehiele that 
ean aehieve LEO using a single stage is seen as a teehnology that ean signifleantly reduee 
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the cost of access to space as well as improving the frequency of access to space. The 
employment of autonomous, reusable launch vehicles requires additional control 
robustness to replace the human’s place in the loop. The guidance and control of these 
autonomous reusable launch vehicles mandate new guidance strategies that are able to 
identify and adapt to vehicle failures during the trajectory and still return to earth safely. 

An important class of failures that concern reentry trajectories are those control 
surface failures which affect the vehicle’s aerodynamic characteristics and controllability. 
The nominal reentry trajectories are designed assuming the vehicle’s original 
aerodynamic characteristics and can no longer be assumed valid following a control 
surface failure. The challenge is to reconfigure the inner loop control algorithm, 
determine the new aerodynamic characteristics of the failed vehicle and to generate a new 
trajectory using the new aerodynamic data to direct the vehicle to a landing point using a 
modified performance criterion. 

This thesis demonstrates the application of a psuedospectral method to determine 
the optimal entry trajectory for the X-33 for nominal conditions as well as for off- 
nominal conditions caused by vehicle control surface failures. Actual X-33 aerodynamic 
data is used in a reduced-order dynamical model to generate the trajectories. Feasibilty, 
optimality and various performance metrics of the psuedospectral approach are 
investigated and quantified. 


B, THE X-33 PROGRAM 

The X-33 was a half-scale technology demonstrator of a potential future Reusable 
Launch Vehicle (RLV) [ref 24]. The prime contractor for the 1.2 billion dollar program 
was Lockheed Martin who called the full-scale vehicle VentureStar. The X-33 was to 
provide data that would lead to the full development of a single stage to orbit (SSTO) 
vehicle that would dramatically lower the cost of access to space through SSTO and 
minimal time between missions. The X-33 represented cutting edge technology in all 
design aspects, from its linear aerospike engines to it composite fuel tanks. The vehicle 
can reach altitudes of 50 miles and hypersonic speeds. The design turn around time for 

this vehicle was 7 days with an eventual turn around goal of two days. The program was 
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expected to begin full flight demonstration in 1999 but was discontinued in 2001 after 
technology maturation difficulties. The X-33 remains a desirable vehicle to study 
because of the large amount of performance data available and because of its multiple 
redundant control surfaces allows for the investigation of a considerable variety of 
recoverable control surface failures and their impacts upon entry trajectories. 

1, X-33 Mission Description 

The X-33 was designed to demonstrate the ability to carry a small payload to low 
earth orbit, return safely to earth, and undergo minimal preparations in order to be ready 
for another mission. The vehicle was to be launched from Edwards Air Force Base and is 
potentially recoverable at a number of suitable fields. The proposed flight test program 
utilized several ascent trajectory profiles terminating at Malstrom Air Force Base, 
Montana or Michael Army Air Field, Fltah. 

The typical X-33 mission begins with a powered ascent that achieves a given 
Mach and altitude at main engine cutoff (MECO). The vehicle performs some activity on 
orbit, transitions and enters the atmosphere passing through several intermediate flight 
stages ultimately ending at landing and rollout. Figure 1 details the various X-33 
mission phases as well as the effectors that are used to control the vehicle. While Figure 
1 shows a reaction control system (RCS) available during the entry phase, this thesis 
assumes that no RCS present during the entry portion of the mission. The weight on 
wheels (WOW) and rollout mission phase are not considered in the trajectory generation. 
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Figure 1 X-33 Flight Phases (from Ref. 5) 


Figure 2 depicts typical X-33 mission profiles with recoveries at both fields. 
Notice that the ascent profiles and the entry initial conditions differ significantly for each 
mission. 
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Figure 2 Representative X-33 Mission Profiles (from Ref. 5) 

A partial list of mission profiles is shown in Table 1. The Maeh number and 
maximum altitude are used as initial conditions for the trajectory optimization routines. 


Flight 

No. 

Trajectory 

Designation 

Landing 

Site 

Dist 

nm 

Mach 

No. 

Max Alt. 
k Feet 

Time 

MECO 

asmg 

Q 

Q 

Alpha 

Fleat 

Rate 

Test Objectives 

1 

9a-2 

Benign 

Michael 

391 

8.84 

186 

239 

4.11 

196 

4864 

1727 

Aero plume interactions. Structural 
interachons, Thrust Vector Control 
Verification Installed Isp 

2 

9d-l 

Transition 

Michael 

391 

10.66 

186 

197 

4.11 

404 

6286 

24.14 

Shock-Shock interactions 

Boundary layer transition 

3 

7b-1 

Max Heatinq 

Michael 

391 

8.97 

189 

221 

4.11 

239 

6576 

27.29 

TPS panel thermal properties 

Aircraft-Like Operational Efficiency 

4 

9c-2 

Real Gas 

Michael 

391 

1021 

189 

204 

4.11 

399 

6006 

22.12 

Real Gas Effects Verification 

Aircraft-Like Operational Efficiency 
Verification 48 hour turn-around 

5* 

7b-l 

Max Heating 

Michael 

391 

8.97 

189 

221 

4.11 

239 

6576 

27.29 

TPS panel thermal propertes / support 
structure 

6 

76 

Real Gas 

Malmstrom 

826 

12.64 

302 

208 

4.03 

342 

7080 

37.05 

Real Gas Effects Verification 
- Surface catalysis with atomic oxygen 

7* 

76 

Real Gas 

Malmstrom 

826 

12.64 

302 

208 

4.03 

342 

7080 

37.05 

- Effects on pitch moment 

8-15 

As Required 

Michael / 
Malmstrom 

391/ 

826 








Repeat Wghts as required to aquire 
additional data 


Table 1. Flight Test Trajectory Matrix (from Ref. 5) 


2, X-33 Vehicle Description 

The X-33 is a lifting body with four sets of aerosurfaces: rudders, body flaps, and 
inboard and outboard elevons. The vehicle also has a reaction control system (RCS) for 
use in unpowered flight at low dynamic pressure. For powered flight, the vehicle uses 
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aerospike engines with thrust vectoring for pitch and roll control and differential thrust 
for yaw control [ref 28]. Figure 3 depicts the general X-33 configuration and identifies 
the vehicle’s control surfaces. 



Figure 3 X-33 Vehicle Depiction (from Ref 1) 

The vehicle gross liftoff weight is 280,000 pounds. The vehicle weight after 
MECO is approximately 78,000 pounds. 

The control surfaces deflection vector, also known as the effector vector, is 
identified by the Greek letter 5 whose elements are identified as 

J = [Relevon,^ Relevon^^^ Lelevon.^ Lelevon^^^ Rudder^j Rudder^^]^ (1-1) 

The X-33 control surfaces have maximum and minimum deflections in degrees as 
shown by the vectors below: 

J^,„=[-30 -30 -30 -30 -15 -15 -60 -30f (1.2) 

^max=[30 30 30 30 26 26 30 60f (1.3) 
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Locked control surface failures are modeled by replacing elements of (1.2) and 
(1.3) that correspond to the failed surfaces, with the positions of the failed surfaces. 

C. IMPLICATIONS OF CONTROL SURFACE FAILURES 

Many published works deal with the subject of reconfigurable control following a 
control surface failure of the vehicle at some point in its mission. The significance of 
reconfigurable control is that it allows the inner loop to compensate for the changed 
aerodynamics of the vehicle and maintain controllability. For many cases, reconfigurable 
inner loop control is insufficient to continue safe flight and a new trajectory based on 
knowledge of the type of failure and the resulting changes in the constraints on the state, 
cost, control and path arguments to the optimal trajectory problem must be computed. 
The replanned trajectory has the ability to continue the mission from the point of failure 
and on to landing thus saving the vehicle and crew which would have otherwise been lost 
as a result of the failure. The availability of an online algorithm such as the AFRL trim 
deficiency algorithm provides the outer loop with the feasible range of Mach and angle of 
attack combinations throughout the trajectory as well as the full 6-DOF effects of the X- 
33 aerodynamic model for use in coefficients of lift and drag calculations for the 
remainder of the trajectory. 

The thesis model uses a pseudospectral method to determine the optimal entry 
trajectory subject to appropriate constraints such as normal load and trim deficiency. The 
approach assumes online identification of the failure. Following the failure identification 
and characterization, the new effector vectors containing the upper and lower control 
surface bounds are passed to the trim deficiency algorithm which determines the new 
feasible ranges of Mach number and angle of attack as well as optimized effector 
displacement vectors for given Mach-dependent angle of attack requests. These values 
are passed to the trajectory algorithm as a path constraint to be used while determining 
the optimal trajectory. An additional function of the trim deficiency algorithm queries 
the vehicle aerodynamic database to provide new vehicle coefficients of lift and drag in 
Mach-alpha space. The outer loop dynamical model employs reduced-order equations of 
motion that use interpolated table lookups to determine the appropriate lift and drag 
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coefficients derived from the vehiele state (based on angle of attaek and Maeh number). 
Onee a eontrol surface failure is deteeted and identified, the aerodynamie table is changed 
to reflect the new vehiele aerodynamies from the perspeetive of the redueed order 
dynamical model. 
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II. APPROACH 


A. GENERAL APPROACH DISCUSSION 

The model used in this thesis uses a reduced-order dynamical model that queries 
the vehicle full order aerodynamic model in order to capture the full 6-DOF model 
characteristics. A flat stationary earth in Cartesian coordinates is used in order to better 
understand and visualize the interactions between the different components of the 
optimization routine. 

B, REDUCED-ORDER DYNAMICAL MODEL 

1. Coordinate Systems 

The simulation uses a three-dimensional Cartesian local horizon coordinate 
system that chooses as its origin the vehicle's position at the beginning of its unpowered 
descent or at the point in the descent at which a failure is identified. The positional 
elements are expressed as pure states in the easily recognized forms of x (downrange), y 
(crossrange), and z (altitude). The coordinate sign convention is the usual xyz sign 
convention where altitude decreasing moving "down", and down range increases to the 
"right" when viewing the xz plane from the side. The y axis is the cross product of the z 
and X axes. Full treatment of dynamical formulations in other coordinate systems can be 
found in Rea [ref. 7]. Figure 4 shows the coordinate system used in the reentry 
simulation. 

where / = flight path angle, the angle made by the velocity vector and the x-y 

plane. 

P = azimuth angle, the angle made by the projection of the velocity 
vector onto the Newtonian x-z plane and the Newtonian x-axis. (This is used 
interchangeably with heading angle since we are dealing with a fiat earth) 

(j) = bank angle, the angle made by the lift vector from vertical. 

V = velocity magnitude 

a = vehicle angle of attack 
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The body axes are depicted in blue; the Newtonian Frame is in black. The 
velocity vector and its projections are red. The lift and drag vectors are colored orange. 
The vehicle is assumed to be in coordinated flight (trimmed in yaw, no sideslip). 

z 



2, Equations of Motion 

The equations of motion are a simplified version of the traditional spherical earth 
equations of motion. There are six equations representing the kinematics and the 3DOF 
dynamics of the problem as adapted from Wiesel [ref. 6]. These equations of motion are 
shown below; 

a. Kinematics: 


x = V cos P cosy 

(1.4) 

y = V sin cosy 

(1.5) 

z = F sin y 

(1.6) 

Dynamics: 


D 

(1.7) 

V =-g sm y 

m 
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( 1 . 8 ) 


Y = 
0 = 


Lcos^ 

mV 

Lsmcj) 

mV 


geos/ 

V 


(1.9) 


The expressions for the lift and drag forces eome directly from 
fundamental aerodynamics, repeated here for completeness. 

i = (1.10) 

( 1 . 11 ) 

In a full 6-DOF model the eoeffieients of lift and drag, 
C^[a,M,(7,S),Cjj[a,M,a,S)avQ functions of several variables: a is the vehiele’s 

angle of attaek, M is the Maeh number, <j is the sideslip angle and 5 is the control 
surface defleetion vector previously discussed. The reduced order model employed by 
this thesis assumes that the coefficients of lift and drag are functions of the state variables 
only. 

C^,C,=f{M,a) (1.12) 

Additionally Mach number and atmospheric density are functions of 

altitude: 


M=M{z) (1.13) 

p = p{z) (1.14) 


This simplified assumption is valid beeause the model performs a table 
lookup of the lift and drag coefficients that were obtained from a full order model that 
uses the vehicle flight condition and the optimized effector displacement vector as 
arguments to interpolate an aerodynamic data table for the vehicle. The resulting 
eoeffieients of lift and drag from wing, body and trim effeets are summed and 
ineorporated into the table. The effeetor displaeement is an implicit argument in the 
algorithm’s eoeffieients of lift and drag determination. The use of the optimized effector 

vector effectively deeouples the outer loop from the inner loop. 
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3, The State Vector 

The state veetor is ehosen to fully represent the state of the system in question at 
any instant. Using the flat earth eoordinate system, the position of the vehiele ean be 
represented by 


X 

T 

z 

where 

X := downrange distanee in feet 

y := erossrange distanee in feet 

z := altitude in feet 


(1.15) 


The veloeity of the vehiele is represented in spherieal eoordinates: 

V := the magnitude of the veloeity in ft/see 

Y := the flight path angle (FPA) in radians 

P := the azimuth angle, in radians 


These three state variables when coupled with the positional information comprise 
the state vector: 


X 

T 


X = 


z 

V 

r 


P 


(1.16) 


4, The Control Vector 

The control vector for this problem is the vector whose elements are the 
controllable parameters of the system. For a vehicle trimmed in yaw, the coefficients of 
lift and drag and therefore x and z motion, are determined by the vehicle's angle of attack. 


12 





The vehicle's bank angle when combined with angle of attack determines the vehicle's 
cross range motion. Since reentry is accomplished dead-stick there is no thrust control 
variable. A suitable candidate for the control vector is therefore 



(1.17) 


a := the vehicle angle of attack 

<j := the vehicle bank angle 

This type of control is also called inertialess control [ref 20] 


C. 4-DOF MODEL 

The inertialess controller can be adapted to more accurately model the physical 
system by recognizing that physical actuators are rate limited. Inertia can be added to the 
system by imposing realistic rate constraints on the actuators. This is accomplished by 
incorporating the elements of the control vector in Equation (1.17) into the state vector in 
Equation(1.16). The new state vector becomes an 8-element vector. 

X 

T 
z 
V 

(1.18) 

r 

P 

a 

<j 

The new control vector components are now the rates of the previous control 
vector’s components in Equation(1.17). 




(1.19) 
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This alternative formulation is provides a more realistie model by preventing 
unlimited eontrol surfaee defleetion rates. Note that while the physieal eontrols remain 
a and cr, the algorithm uses the rates of these variables in its optimization. This 
dynamieal model is ealled a “4-DOF” model beeause it provides additional fidelity 
beyond the 3-DOF model. The full 4-DOF trajeetory optimization algorithm is shown in 
Figure 6 


It 


li 



Error 


Optimization Routine 


Figure 5 Trajeetory Optimization Feasibility Error Estimation 

The magnitude of the error shown in Eigure 5 is used to ehoose the size of the nodes 
eomputed in the psuedospeetral method. The number of nodes is inereased until the error 
between the propagated states and optimal states is aeeeptable. Aeeuraey eomes at a 
eost; the larger number of nodes signifieantly inereases algorithm exeeution time. 


D, ATMOSPHERIC MODEL 

An AERE atmospherie model whieh generates speed of sound and air density for 
a given an altitude is used. The outputs of the atmospherie model are eombined with the 
vehiele veloeity to determine the vehiele’s Maeh number to be used in the aerodynamie 
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table lookups as well as the control allocator and trajectory path constraints. Figure 6 
and Figure 7 depict the model’s variance of the speed of sound and atmospheric density 
as a function of altitude. The atmospheric model is a combination of exponential 
functions separated by discrete altitude breakpoints. 
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Figure 6 Plot of Density vs. Altitude (from AFRL atmosphere model) 
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Figure 7 Plot of Speed of Sound vs. Altitude (from AFRL atmosphere model) 

E. GRAVITATIONAL MODEL 

The gravitational model used is a low order model that uses a constant 
gravitational acceleration which is assumed to always act in the negative z direction. The 
decision to use a simplistic model is based primarily on the desire to remove additional 
complications and observe the simulation dynamics in as pure a manner as possible. 


F. X-33 AERODYNAMIC MODEL 

This thesis considers the atmospheric entry of the X-33 reusable launch vehicle 
following a mission in low-earth orbit. The computed trajectory is an unpowered descent 
to a known landing field. The trajectory exposes the vehicle to full ranges of Mach and 
angle of attack. To accurately model the aerodynamic performance of the X-33, the 
trajectory code utilizes a dynamic linked library (dll) developed by Dr. David Doman of 
the Air Force Research Laboratory (AFRL) consisting of compiled C++ code that queries 
X-33 aerodynamic data derived from the full MAVERIC simulation. The file, 
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x33aero.dll, determines the full order aerodynamie eharaeteristies ineluding wing-body 
and trim effeets from an input veetor eonsisting of Maeh number, angle of attaek and 
eontrol effeetor position. Sinee the position of the eontrol effeetors is arbitrary for any 
given desired angle of attaek, two-dimensional tables of total lift and drag eoeffieients in 
Maeh-alpha spaee are generated using another AFRL eode that uses a pieeewise linear 
programming (PLP) algorithm to determine the optimal eontrol alloeation of the vehiele's 
eontrol surfaees (i.e., eontrol indueed moments balanee the moments generated by the 
wing-body with minimum eontrol defleetion in a 1-norm sense) to aehieve a desired 
angle of attaek. The optimal eontrol alloeation is not direetly used by the outer loop 
trajeetory generator, but the resulting eoeffieients of lift and drag derived from the 
optimized eontrol effeetors for a given eommanded angle of attaek and Maeh are used in 
the outer loop. 

1. Mach Breakpoints 

The X-33 aerodynamie data eovers a range of angles of attaek in inerements of 
one degree from -10 to 50 degrees over a diserete Maeh range given by [0 0.3 0.6 0.8 0.9 
0.95 1.05 1.15 1.46 2.01 2.75 3.49 4.75 6 10 20]. This diserete Maeh interval is used in 
the eompiled C++ eode refereneed by the DIDO dynamies funetion. The aerodynamie 
database output is in the body frame, so an appropriate transformation to a Newtonian 
frame is neeessary in order to analyze the trajeetory. The normal foree path eonstraints 
ean be analyzed in the body frame. 

2. AFRL PLP Control Allocation Model 

The AFRL eontrol alloeation algorithm is used to trim the vehiele as well as 
eomputed the lift and drag estimates for the lift in any eonfiguration. The algorithm 
assumes that the vehiele is in symmetrie flight (sideslip angle is zero) and the resulting 
lateral wing-body foree and moment eoeffieients are zero. By assuming a zero magnitude 
body-axis angular veloeity, the wing-body pitehing moment of the vehiele beeomes a 
funetion of Maeh number and angle of attaek. This allows for the neeessary eonditions to 
rotationally trim the vehiele that the moments resulting from the eontrol surfaee 
defieetions must be equal and opposite that of those produeed by the wing-body [ref. 22] 
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= 
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Where l,m,n denote the roll, piteh and yaw axis, is the moment coefficient 

caused by the control surface deflection about either the l,m ox n axis. M is the Mach 
number, a is the angle of attack and ^ is the control surface deflection vector caused by 
the eight individual vehicle control surfaces. 

The first task of the control allocator determines the optimum effector 
displacement that results in minimum control deficiency. The full description is 
contained in reference 22 and repeated here for completeness. The full set of control 
surface deflections for the X-33 can be represented by the eight-element effector 
displacement vector previously identified in equation (1.1) . This vector is subject to the 
limits of equations (1.2) and (1.3) or a subset thereof for control surface failures. The 
control allocation problem uses an iterative algorithm that queries the full 6-DOF X-33 
aerodynamic database. The piecewise linear constrained control allocation problem is 


MirigJ^ = 


[M,a,d) 4- {M,a,d) + C„^, {M,a,d) + (M,«) 


subject to 


Notice that the upper and lower bounds on d are dependant upon the operating 
mode of the vehicle. A vehicle operating under nominal conditions has the full range of 
deflections available while a vehicle operating with a control surface failure has a subset 
of the nominal displacement vector’s range. 

A nonzero value of means there is insufficient control power to rotationally 

trim the vehicle (it is trim deficient), and the vehicle will depart controlled flight. The 
trim deficient regions in Mach-alpha space can be mapped and identified as “keep-out” 

zones for the outer loop trajectory algorithm. If the control allocation algorithm is able to 
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rotationally trim the vehicle, any excess control power can be used to maneuver the 
vehicle. The control allocator employs a second control allocation optimization 
algorithm to minimize the deviation from a minimum trim-drag condition. This 
condition is represented by an eight-element preference vector set to zero. The second 
optimization is 


MingJ^ 




subject to ^ ^ ^ 

C,AM,a,sU<i 


In the case where the vehicle is operating with excess control power, the effector 
displacement vector is the output of the second optimization routine, otherwise ^ 
minimizes the trim deficiency. The effector displacement vector that is the solution of 
the control allocator is represented as ^ to delineate it from the infinite combination of 
control surface deflections that are available to achieve the same desired control but are 
not optimized for either trim deficiency or minimum trim-drag. 

Given S* [a, M) (from the control allocation algorithm), the X-33 aerodynamic 

database is queried over the Mach-alpha grid to generate the wing-body and trim lift and 
drag coefficients. These coefficients when summed together represent the total lift and 
drag coefficients. The use of the control allocator solution effectively decouples the 
trajectory generation algorithm from the inner loop and allows the lift and drag 
coefficients to be functions of angle of attack and Mach number. Effector displacement 
remains an implicit argument and allows for the full 6-DOF model’s behavior to be 
captured by a lower order model. The outer loop equation (1.12) is elaborated using the 
lift and drag coefficients as the sum of the wing-body coefficients plus the effector- 
dependent trim coefficients using the control allocator solution. 
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Q(M,a) = Q (M,a) + Q^,(M,«,^’ {M,a)) 
Co{M,a) = Co{M,a) + Cj,^(M,a,S* {M,a)) 


( 1 . 20 ) 


The relationships in (1.20) hold for a vehiele in symmetric flight and can be 
stored in an aerodynamic lookup table for each Mach and angle of attack combination 
subject to the limits on the effector displacement vector dictated by the presence of 
nominal or failure conditions. Each failure condition requires a separate run of the 
control allocator. A block model implementation of the control allocator is depicted in 
Figure 8 



Aerodata for 


^oplinid 


Figure 8 AFRF Control Allocator 

3, Trajectory Replanning Algorithm 

The vehicle is assumed to operate in the nominal condition until a failure has been 
identified through some online fault detection or identification algorithm at which point 
the lift and drag table corresponding to the failed vehicle configuration is generated. An 
approximate block diagram of the guidance loop is shown in Figure 9 


20 






Figure 9 Trajectory Replanning 

4, Failure Mode 

Once a failure has been identified, the PLP control allocator generates a new set 
of control deflections for the failed vehicle and uses the aerodynamic database to 
determine the new trim lift and drag contributions. In the process of computing the 
control deflections for the failed vehicles, points in the Mach-alpha envelope for which 
the vehicle cannot be trimmed are identified. These trim deficiency regions must be 
avoided and are thus used as path constraints for the trajectory replanning. The trajectory 
aerodynamic data table is changed to one appropriate for the vehicle in the current failed 
condition and a new trajectory is generated for the new aerodynamics and path 
constraints. 

5, Failure Mode Effects 

The trim deficiency path constraint on the available range of angle of attacks for a 
given Mach number requires that the vehicle must be rotationally trimmable (assuming 
no side slip) for a given angle of attack. When the vehicle operates in a nominal 
configuration it is capable of trimming over a full range of angle of attack, however, 
when the vehicle experiences control degradation, its available range of angles of attack 
may be severely reduced. Doman, et al. [ref 19] give full treatment of this subject. Figure 
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10 depicts the available values of angle of attack that will allow for a rotationally 
trimmed vehicle experiencing no failures. Notice that because the maximum trim 
deficiency value is on the order of machine precision so that the entire range of 
combinations of angles of attack and Mach numbers is considered to be available in the 
nominal (unfailed) case. 
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Figure 10 X-33 Trimmable Values for Angle of Attack.(from Ref 19. ) 

Figure 11 shows an example of the constrained Mach-alpha envelope available to 
the X-33 following a dual body flap failure (stuck at 26 degrees up). Notice that the 
envelope describing allowable angles of attack (in dark blue) previously depicted in 
Figure 10 is now severely reduced. Angles of attack other than those contained within 
the dark blue region will place the vehicle in an attitude where it is not trimmable in pitch 
(pitch deficient) and therefore uncontrollable. This establishes a restricted region that 
feasible trajectories must avoid. In other words, the trajectory must have a clear path in 
Mach-alpha space from the initial conditions to the final conditions. A clear path is 
defined as a path from one point defined by the flight condition (Mach-alpha) to another 
point (Mach-alpha) with no intervening regions of unacceptable trim deficiency. Should 
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a path cross into or through a region of pitch deficiency, the trajectory becomes infeasible 
because the vehicle is no longer rotationally trimmable and therefore uncontrollable. 


Pitch Deficiency Contour 
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Figure 11 X-33 Trimmable Values of Angle of Attack following a Dual Stuck Flap 

Failure (failed at 26 degrees up), (from Ref 19). 

A feasible solution, if it exists, must be found using these new path constraints. 
The allowable range of angle-of-attack can be determined in a manner similar to the 
previous method used to generate the lift and drag tables using the AFRL X-33 
aerodynamic database and control allocator. The result is a dataset that accurately 
describes the trim deficiency as a function of Mach number and angle of attack (the 
vehicle’s flight condition). This dataset is the output of the AFRL trim deficiency 
algorithm which also determines the new coefficients of lift and drag for the failed 
condition using the new optimized effector displacement vector. 

Some failures will have very restrictive path limits for feasible trajectories while 
others may have no clear path. Figure 12 shows a failure where the vehicle is free to 
operate in a high Mach flight condition but cannot safely land. Figure 13 depicts a 
failure where the vehicle is limited to slower speeds and alpha ranges to safely land. 
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Note that the previous diseussion does not consider any other limitations such as normal 
force, heating or dynamic pressure on the trajectory. These constraints will further 
reduce the area available to a clear path for a feasible trajectory. In some cases, the 
additional constraints will completely eliminate a clear path and render trajectory 
infeasible. It is easier to visualize the optimal trajectory with constraints as a region 
through Mach-alpha space for which the vehicle must be rotationally trimmable. If a 
failure impinges on this corridor and prevents a clear path through to landing, then there 
are no feasible trajectories available. 
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Figure 12 X-33 Pitch Deficiency Following a Stuck Right Rudder (15 degrees) 
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Figure 13 X-33 Pitch Deficiency Following a Stuck Flap (at -15 degrees) 

6, Modeling Trim Deficiency 

As the previous graphs of failure modes show, the trim deficiency can be modeled 
as a function of two variables; 

d = f{M,a) (1.21) 

The function, f, can be determined from various fits to the trim deficiency data. 
This thesis uses the Table Curve 3D software, version 4.0 by Systat Software 
{www.systat.com) to generate MATLAB function codes that approximate the data 
generated by the trim deficiency routines. Table Curve 3D is very robust and can apply 
and analyze the fits of over 9000 equation types to the fitted data in short order. An 
example of the raw data from the trim deficiency code imported into Table Curve 3D is 
shown in Figure 14 
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Trim Deficiency in Mach-Alpha Space 
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Figure 14 Table Curve 3D Trim Deficiency Data for a Dual Stuck Flap Failure (26 

degrees) 


A 10* order Chebyshev polynomial produces a very good data fit. The function 
generated by the Table Curve 3D code was used in the path function tile of the DIDO 
code. The path function file combined with appropriate bounds on the path prevents the 
selection of combinations of Mach and angle of attack that exceed a given trim 
deficiency. An important selection criterion for the fitted curve is the function’s behavior 
at the boundaries of the Mach-alpha envelope. Some fitting functions that showed good 
residuals for the data were rejected out of hand because the low Mach and angle of attack 
data tits were very poor. It is more important that the fitting function models the data 
within the feasible ranges of trim deficiency than for ranges where feasible trajectories do 
not exist. Mismatched data for extreme trim deficient regions can be tolerated but it is 
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essential that the fit near the transition from zero trim deficiency to positive trim 
deficiency be as good as possible. 


Trim Deficiency in Mach-Alpha Space 
Rank 1 Eqn 449 Chebyshev X,LnY Bivariate Polynomial Order 10 
r2=0.99796074 DF Adj r2=0.99790679 FitStdErr=0.0008263599 Fstat= 18791.943 
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Figure 15 Trim Deficiency Fitted with 10th order Chebyshev Polynomial 

The trim deficiency can also be modeled as a pure table lookup which queries the 
tablular data from the trim deficiency algorithm. This approach is more accurate than the 
fitted function but was omitted from this study because of the inherent execution 
penalties resulting from a 3-dimensional table lookup. The feasibility of following a trim 
deficiency path constraint is shown by projecting the trajectory into Mach-alpha space 
along with the trim-deficiency contour plot and verifying that the trajectory does not 
violate the path constraints. The feasibility can also be shown by inspection of the 
numerical trim deficiency data for the trajectory. 
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III. THE OPTIMAL CONTROL PROBLEM 


A. OPTIMAL CONTROL THEORY 

The following sections summarize the general optimal control problem 
formulation and the two known methods of generating solutions. 

1. Preliminaries 

Nearly any trajectory generation problem can be posed as an optimal control 
problem and the problem considered by this thesis is no exception. Given a system with 
dynamic constraints 

i = /(x,M,r) (1.22) 

where 

/ := vector of functions that describe the dynamics of the system 

X := vector of states that fully describe the system at any r 

u := vectors of control variables 

T := an independent variable (usually time) 


The system is subjected to additional constraints to include the path constraint: 
hj < h[x,u,t)< h^, (1-23) 

where 

h := vector of functions describing the path constraints 

hj := vector of lower path bounds 

:= vector of upper path bounds 
and boundary conditions 

5 ^e(x(ro),x(r^),ro,r^)<^ (1.24) 


where 
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vector of initial boundary conditions (at ) 


e{x{T^),T^) := 

:= vector of final boundary conditions (at ) 

e^ := vector of lower bounds 

:= vector of upper bounds 

The system also has bounds on the control and state variables represented by 



( 1 . 25 ) 

^ <M(r)<^ 

( 1 . 26 ) 


where 


x{t) := 

state vector at any instance of r 

M(r) := 

control vector at any instance of r 

:= 

vectors of lower and upper state bounds 

:= 

vectors of lower and upper control bounds 


The bounds on the state and control vectors can be thought of connecting the 
physical system being modeled to the mathematical model used to solve the optimal 
control problem. For example, if the control is acceleration, no physically realizable 
device produces infinite acceleration and we therefore may choose to limit the control to 
a sensible value. Similarly, if one element of the state vector is velocity, we may choose 
to limit the velocity to one that is within the design limits of the physical system. 
Imposing constraints on the state and control vectors may not be necessary in all cases 
but the resulting optimal solution of unbounded states or controls needs to be given a 
hard reality check when realizing a physical implementation of the optimal solution. 

An optimal control problem seeks to determine the solution of the preceding 
dynamical system with its bounds and constraints to minimize a given performance 
index. The performance index is mathematically defined below: 
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where 


y(x(r),M(r),ro,r^) = £'(x(rj,x(?^.),r„,r^) +j7^(x(r),M(r),r)(ir (1.27) 

^0 

E := scalar cost function evaluated at the boundaries 

F := scalar cost function evaluated over the entire interval 

When J consists of only E, the cost function is in Mayer form. When J consists of only 
F, the cost function is in Lagrange form. When both E and F are present, the cost 
function is in Bolza form. 

The general optimal control problem is fully posed in the following manner; 

min„ 7(x(r),M(r),r„,r^) = £'(x(rj,x(t^),r„,r^) +jF(x(r),M(r),r)(ir 
subject to x = f (x, u, r) 
hi - h[x,u,T) < 

5 ^e(x(ro),x(r^),ro,r^)<^ 

Ui<u{t)<u^ 

The formulation of the system as represented in the state vector, the appropriate choice of 
the dynamical equations, and the selections of the constraints and bounds forms the basis 
of the fully posed optimal control problem 

2. Solution Methodology 

Once the problem has been properly formulated there are two general methods 
available to solve the optimal control problem, direct and indirect [ref 27]. Indirect 
methods tend to generate fast solution times and greater accuracy but are considerable 
more difficult to formulate and are very sensitive to the initial guess [ref 27]. 
Formulation of the indirect method can be exceedingly difficult in situations where the 
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dynamics functions are not pure functions but instead rely on numerical data such as the 
aerodynamie lookup tables used in this thesis. 

Direct methods reduce the optimal control problem to a single large Nonlinear 
Programming (NLP) problem. The strengths of the indirect methods are that the 
formulation is significantly easier and the methods are relatively insensitive to the initial 
guess. A partieular drawback of direct methods is large exeeution times. When direct 
methods are used, there are tradeoffs between exeeution time and aecuraey. An inerease 
in the size of the NLP program does not necessarily yield a proportional inerease in 
solution accuracy or solution time. This drawback inhibits real time optimal control using 
direct methods for certain problems such as are considered in this thesis; however, as 
computing power continues to increase, the potential for real-time optimal control will 
soon be realized. 

3. DIDO 

This thesis uses the DIDO numerieal dynamie optimization software developed 
by I.M. Ross and F. Fahroo of the Naval Postgraduate School. DIDO employs a direct 
Legendre psuedospectral technique that uses the NLP solver SNOPT. The reader is 
referred to References 12 through 17 for a full description of the Legendre psuedospectral 
method. This thesis shows another application of DIDO in solving optimal trajectory 
problems. 

4, Optimality 

Since the formulation of a problem using a direct Legendre psuedospectral 
method does not determine the adjoint equations, another tool must be employed to link 
the solutions generated by indirect methods and direct methods. The Coveetor Mapping 
Theorem (CMT) provides this link and is fully explained in Reference 8. The 
implications of the CMT allow a determination of the optimality of a given solution by 
comparing the costates of the solution to the Kamsh-Kuhn-Tucker (KKT) multipliers. 
Every time a state or control becomes bounded, the corresponding costate should follow 
the KKT conditions as if they were the KKT multipliers. Additional evidence of 
optimality can be obtained by investigating the solution's Hamiltonian and the final value 


32 



of the cost function. These necessary conditions may be easily obtained from the DIDO 
solution and provide insight into the optimality of the DIDO generated solution. 

5, Feasibility 

Because the DIDO solution is based upon the controls for a discretized set of 
points determined by the Legendre polynomials used in the solution, the DIDO solution 
may not possess the desired accuracy at the final endpoint. Since the accuracy of a DIDO 
solution is a function of the number of nodes evaluated in the solution space, all DIDO 
solutions should be evaluated for feasibility. This is usually accomplished by comparing 
the DIDO solution to a separate solution using the DIDO controls propagated by a Runge 
Kutta ODE solver such as MATLAB's ODE45. The error between the propagated states 
and the DIDO-derived states gives some measure of confidence of the accuracy of the 
DIDO solution. The required accuracy is determined by the application, for example 
while the downrange accuracy of a 20-node DIDO solution reentry trajectory is 0.0001% 
of the propagated downrange distance, it still may be inaccurate enough to prevent a 
feasible landing at the desired field. Eeedback may be able to overcome this inaccuracy 
provided the solution update is achieved at a reasonable rate. Another guard against 
inaccuracies in a trajectory is to provide additional margins when defining the endpoint 
conditions. 

6. Scaling and Balancing 

While direct methods are relatively insensitive to the initial guess, a well- 
formulated problem can be completely ruined by improper scaling and balancing of the 
variables and constants used in the computation. Ross [ref 21] notes the effect of 
improper scaling on the DIDO solution. It should be noted that inks units are a much 
better choice of units and are preferred to English units from a computational viewpoint 
because mks units are direct combinations of distance, mass and time base units and thus 
easily scaled while some English units have no direct equivalent. 

The intent of scaling and balancing is to avoid the problems associated with 
variables whose range exceeds that of numerical precision (about 16 significant figures 
on a Pentium-class PC). A proper scaled problem scales the variables to a range of 0 to 
1. While it is not necessary to be exactly within this range, it is clear that variables of the 
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same order of magnitude avoid problems assoeiated with eomputing ill-eonditioned 
matriees resulting from numerieal impreeision. The general approaeh to sealing the 
problem variables is as follows: 

The fundamental units in any system of measure traditionally eover mass, time 
and distanee. Other units sueh as area, volume and veloeity are derived from these 
fundamental units. Consider the initial eondition for the third element of the state veetor, 
altitude. The initial altitude is 125,000 feet while the final altitude is 500 feet. 
Additionally, the initial downrange distanee is at the origin of the eoordinate system 
while the expeeted maximum downrange is nearly 4,000,000 feet. A sealed distanee unit 
is ereated by applying a scale factor to the state variables. 


X = 


X 

y 

DU 


z = 


z 

DU 


Where DU represents the sealing of the distanee unit. To seale a variable that 
represents veloeity notiee that 


, . distanee 

velocity =- 

time 

To seale veloeity, simply divide by the veloeity sealing faetor. 


V = 


V 

VelUnit 


V 

'DU^ 

[YuJ 


= v 


' TU'" 
^~DU, 


Similarly, the sealing faetors for aeeeleration and foree are 


AcelUnits = 


DU 
{TUf ■ 


ForceUnits = MU 


DU 

Stu)" 




J 
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Arbitrarily choosing DU=100,000 ft and TU= 50 seconds for this problem seales 
the state variables as shown: 


X. 


yf 


z 


f 




0 


Pa- 


“ 0 “ 

0 


A 


0 

200,000/t 


f. 


1 

8000 ft! see 




2 


4,000,000/r 


A 


“20 

0 


A 


0 

500/t 


A 


0.0025 

250/t/see 




0.0625 


After sealing, the range of distanees falls between 0 and twenty for down range 
and between 0 and 2 for veloeity. The problem should also be balaneed, that is all path 
eonstraints should be as elose to the same order of magnitude. While this may not be 
possible, every attempt should be made to balanee the problem as well as to seale it. 

B, REENTRY VIEWED AS AN OPTIMAL CONTROL PROBLEM 

After sealing and balaneing, the atmospherie entry (or reentry) ean be formulated 
as an optimal eontrol problem. 

The reentry problem eonsiders a vehiele entering the atmosphere from an initial 
position, presumably near LEO and seeks to generate an optimal trajeetory that brings it 
to a final state speeified by an endpoint eondition. The trajeetory may not violate the 
path eonstraints whieh limit heating rate, normal foree, and axial foree (dynamie) 
pressure. The optimal trajeetory seeks to minimize a performanee index whieh is 
generally a funetion of the state of the vehiele. Typieal performanee indiees for reentry 
are maximizing downrange, maximizing eross range, minimizing eontrol effort, or some 
weighted eombination of these indices. 
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The trajectory is also subject to dynamic constraints which describe the dynamics 
of the system. As discussed in the previous section, this thesis uses a simple dynamic 
model that retains the benefits of a higher fidelity model by including the effects of trim 
lift and drag and using table lookups into a high-fidelity aerodynamic database. 

For a typical trajectory that seeks to determine the maximum downrange value for 
an entering vehicle, the cost function is simply a Mayer cost: 

j(x(t),M(t),to,t/) = -X/ 

Time is the independent variable for this problem. A Lagrange component of the 
cost can be included as desired. The initial boundary conditions for a representative 
mission are 


X. j 


“0 

L- 


0 



125,000/t 



Mach 8 

A 


-1.3 deg 

A 


0 

a, 


0 



0 


The endpoint equation does not need to fully describe the state vector, only the conditions 
that are necessary for a safe transition and landing. The manifold can be simple such as: 


The vector above does not specify the final angles of attack or bank and it is possible that 
the vehicle could arrive at the ending manifold in an unrecoverable position (high angles 
of attack and bank along with large vertical speed). To prevent this from occurring, it is 
good practice to fully specify these variables as well as the associated acceptable range of 
values at the end point. Initial trajectories computed appear to be feasible trajectories but 



) _ 



“500 ft 


_^f\ 


Mach 0.15 
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upon closer inspection met the end point manifold with vertical speeds exceeding 20,000 
feet per minute. The endpoint manifold was then modified to an endpoint set limited to 
vertical speeds between -1500 ft/min < z < 500 ft/min . 

The path constraints are derived primarily from the design and mission considerations of 
the vehicle. 
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The limits that have values of plus and minus infinity represent states without 
hard limits and are not truly constraints upon the state vector. The constraint on the 
altitude component of the state vector is imposed to prevent premature contact with the 
ground. 
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IV. APPLICATIONS 


A. FOOTPRINT DETERMINATION 

The first step in determining the vehiele performanee begins with investigating 
the nominal (unfailed) ease. A rough vehiele footprint ean be established by ehanging 
the eost funetions to inelude various eombinations of the form 

M = Ax(^end) +Byi^end) 

Where A + B = 1 and 0<y4<l;0<5<l. 

x[end) and y[end) are the final values of downrange and erossrange. 

The eonstants for eost funetion used to determine the maximum downrange are 
simply =7 it follows, therefore that B=0. Note that when y4=-7, minimum downrange is 
implied. The eonstants for eost funetion used to determine the maximum erossranges are 
simply B=1 it follows, therefore that A=0. Note that when B=-l, minimum erossrange is 
implied, whieh for the nominal ease should be symmetrie about the x axis. 

B, LANDING AT A SPECIFIC FIELD WITHIN THE FOOTPRINT 

Onee the rough limits of the footprint have been established, the vehiele should be 
able to land at any field within the footprint. In order to reeover the vehiele at a speeifie 
field, the eonditions at the landing threshold that lead to a “sueeessful” landing must be 
strietly defined. Original runs that speeified a small subset of the state veetor at the 
endpoint led to unsatisfaetory eonditions at the landing threshold. Consider the endpoint 
eonditions speeified by 

e(xf,tf) = 

Optimal trajeetories ean be generated whieh meet these endpoint eonditions; 
however, the vehiele state may be sueh that a sueeessful landing is not possible given the 
praetieal limitations of the vehiele. Original trajeetories whieh used the previous 
endpoint eondition appeared reasonable until the vertieal speed at the landing threshold 
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was investigated. Some trajeetories hit the landing threshold with vertieal speeds in 
exeess of 4000 feet/see. Clearly, this is an unreeoverable attitude- the vehiele passes 
through the endpoint on its way to its final resting plaee at the bottom of a very deep 
erater instead of the end of the runway. Another eonsideration is angle of bank. Hitting 
the terminal eonditions with strong left drift and extreme right bank angle ean be 
disastrous. Limits on the bank angle are elearly needed. The endpoint eonditions that 
speeify a landable attitude have now grown to 




Vf 


Notiee that angle of attaek is not speeified in the endpoint eonditions and there is 
a ehanee that the vehiele eould be near stall. The endpoint eonditions as stated above are 
insuffieient to hit a speeified point along the trajeetory. The landing field is speeified by 
Cartesian eoordinates x,y,z as well as a runway heading, speeified by the azimuth angle, 
P . Notiee that in the flat earth model, the altitude of the field is assumed to be zero but 
is not neeessarily zero. Endpoint eonditions for landing at fields with nonzero elevations 
are simply aeeommodated as 

~ ^field ^ojfset 

The new endpoint eonditions for landing at a speeified field in a landable vehiele 
attitude are now speeified as 

e[xf,Xf,tf) = [_Xf z^ 


The runway heading for the seleeted landing field ean have signifieant effeet on 
the trajeetory beeause reversing direetion requires eonsiderable energy loss during the 
trajeetory. It is entirely possible that a footprint generated with loosely speeified 
endpoint eonditions does not eontain interior points for feasible trajeetories when the 
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endpoint conditions are more fully specified (such as a runway heading 180 degree from 
the original vehicle azimuth). 
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V. RESULTS 


A. DISCUSSION 

The following sections discuss the footprint and individual trajectory analysis for 
several cost functions. Unless otherwise specified, the trajectories for both failed and 
nominal vehicles use the same initial conditions of 

e(x,,t,) = [x, y, z. V, y, a, 

With the following values 

e(x.,t.) = [0 ft Oft 125 kft Mach 8 -1.3° 0° 0° 0°f 

These initial conditions correspond to a point on the entry trajectory after the peak 
heating and dynamic pressure and are suitable for failure as the vehicle can begin to use 
its control surfaces to control the trajectory. The initial condition for the angle of attack 
may be too shallow based on previous work, but it is of little consequence for studying 
the effects of failure on trajectory replanning and retargeting. 


The endpoint conditions for both trajectories are of the form 

= Zf Zf (pj-'J 

With the following values 

= [free free free -500±1000 ft/min Mach 0.15±0.003 free 0±20°]^ 

Both vehicles observe a normal force constraint of 3±/- 6 g’s. The vehicle experiencing 
a control surface failure must observe the trim deficiency path constraint while the 
nominal vehicle is unconstrained by trim deficiency because as shown in Figure 10 , the 
vehicle has adequate control power to trim over the entire flight envelope. 

Once the rough footprint boundaries are established for the vehicle experiencing 
control surface failures, the algorithm’s ability to generate trajectories to a specific 
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landing field and its associated runway is demonstrated. For this case, the Cartesian 
coordinates and azimuth are no longer free and are replaced with the field xy coordinates, 
an altitude corresponding to 500 feet above the ground and a corresponding runway 
azimuth. 

B, FOOTPRINT ANALYSIS FOR NOMINAL AND STUCK FLAP 

TRAJECTORIES 

An approximate footprint for the nominal and stuck-flap trajectories can be 
determined by choosing the cost functions corresponding to minimum and maximum 
down range distance as well as minimum and maximum cross range distance. The four 
points form an approximate boundary of the footprint and feasible trajectories to the 
interior points of the polygon formed by these boundary points should exist. Note that 
feasible trajectories to exterior points are not ruled out but should be proven for each 
case. Additional boundary points can be determined by varying the A and B terms in the 
previous equation M = Ax [ end) +By [ end ). 

1, Aerodynamic Considerations 

The introduction of a control surface failure can drastically change the 
aerodynamics of the vehicle. The coefficients of lift and drag for the vehicle 
experiencing a stuck flap have changed significantly. The constraints on angle of attack 
resulting from the trim deficiency requirements produce a vehicle which may look the 
same but flies considerably different. The footprints for the failed and normal trajectories 
should appear different because of the different aerodynamic configurations and Mach- 
alpha combinations. 

Despite the aerodynamic changes of the failed vehicle, its performance should be 
a subset of the nominal vehicle’s performance. This is a valid assumption because the 
control surface failures encountered in flight limit the control surface deflections to a 
subset of their nominal values. In other words, the nominal vehicle can fly the same 
profile as a vehicle with control surface failures simply by voluntary limiting the unfailed 
control surface deflections to those of the vehicle with the failures. The only difference 
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becomes semantic; the nominal vehicle “simulates” the failure while the vehicle with the 
actual failure has no choice but to fly (if it can) with the current failed conditions 

2, Footprint Observations and Conclusions 

Figure 16 shows the footprints for four selected trajectories: nominal, stuck flap 
failure with loosely specified endpoint conditions, stuck flap failure with tightly specified 
endpoint conditions and stuck flap failure with tightly specified endpoint conditions and 
trim deficiency path constraint. Several things are readily apparent: 

a. The footprints for vehicles with failures have better minimum downrange 
performance than that of the nominal vehicle. This was surprising but can be explained 
by the control allocation algorithm. Recall that for a given Mach-alpha pair, the control 
allocator computes a vector of control deflections that simultaneously balance the 
moments produced by the wing-body and minimize the 1-norm of the control deflection 
vector. Changing the preference vector to maximize trim-drag or control deflections 
within the control allocator can solve this particular issue; however, the differences in the 
footprints are indicative of greater problems associated with having the outer loop 
decoupled from the inner loop. 

b. There is great variance between the three footprints of the failed vehicles. 
The largest footprint for the failed vehicles represents the traditional footprint for a failed 
vehicle as found in the literature. One such approach to determining this footprint uses 
Vinh’s [Ref 23] equations. The footprint uses loosely specified end point conditions 
(altitude and airspeed), and does not prevent the angle of attack-Mach combination from 
entering regions of the flight envelope where the vehicle cannot achieve moment 
equilibrium due to a lack of control power. As such, the larger footprint depicts where 
the “pieces fall” since in the early moments following a failure, the trajectory would 
violate the trim deficiency constraint and the vehicle would become uncontrollable and 
break up. 

c. The effect of more restrictive endpoint conditions further reduces the 
footprint. The third largest footprint represents the footprint for the trajectory for the 
failed vehicle that observes the trim deficiency path constraint. While this trajectory 

appears to be sufficient, further inspection reveals that the endpoint conditions are met 
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with a vertical speed of nearly 20,000 feet per minute. The vehicle meets the endpoint 
conditions on its way to a deep crater as there is insufficient remaining altitude to slow 
the rate of descent to an acceptable level. The smallest footprint represents a trajectory 
for the failed vehicle with tightly specified endpoint conditions that limit the vertical 
speed and angle of bank to levels that permit landing. This footprint for a “landable 
vehicle” in the failed condition is 1% of the area of the largest failed footprint and shows 
that trajectories generated using energy analysis with loosely specified endpoint 
conditions and without trim deficiency path constraints result in overly optimistic 
footprint estimates. 


Footprint Comparisons of Various Trajectories 



Figure 16 Comparison of Nominal and Stuck-Flap (+26°) Footprints 
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Footprint Comparisons of Various Trajectories 



Figure 17 Contraction of Failure Stuck-Flap (+26°) Trajectory Footprint 


The difference between a “landable” and “unlandable” footprint is shown in 
Figure 17 Notice that the introduction of vertical speed into the endpoint set reduces the 
footprint noticeably. Landable footprints result from trajectories with endpoint 
conditions that allow a safe transition from the landing threshold to weight on wheels and 
rollout that are within the vehicle and payload performance parameters. 

3, Maximum Downrange 

The following trajectories for the nominal and failed vehicle were conducted as 
part of the rough footprint analysis. The performance criterion for these trajectories was 
to maximize the downrange position at the end of the trajectory. The following trajectory 
analysis is a detailed comparison of the nominal and failed vehicle trajectories for each 
associated cost function. 
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a. Mach-Altitude Profile 

The Mach-altitude profiles for both trajectories are shown in Figure 18 . 
Notice that the nominal trajectory trades excess speed for altitude. At the top of the zoom 
climb, the vehicle assumes the best lift/drag configuration and continues to the surface. 
The trajectory for the vehicle with the stuck flap is prevented from using the same 
combination of angle of attacks as the normal because of the imposed trim deficiency 
requirement. The trajectory does use the largest allowable angle of attack but this is 
insufficient to produce the zoom climb observed in the nominal trajectory. 



Mach number 


Figure 18 Maximum Downrange Mach-Altitude Profile 
b. Trajectory in Cartesian Coordinates 

Figures 19-22 depict the Cartesian coordinates of the trajectory. The x, y, 
and z coordinates are direct results of the controls and the time of their applications. The 
only performance measure for these trajectories is to maximize the downrange (x). There 
is no limit placed on crossrange (y) and the altitude must meet the endpoint condition of 
500 feet. 
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The downranges for the failed and nominal trajeetories differ by nearly 
200 miles. This is not unreasonable considering that the symmetrical flap failure 
drastically increases the drag on the vehicle. This drag when combined with the imposed 
trim deficiency path constraint considerably slows the vehicle and prevents it from 
trading altitude for velocity. The total trajectory time of flight is much shorter. 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 19 Downrange Profile for Maximum Downrange Trajectories 

The crossrange profiles of both trajectories exhibit some unexpected left 
drift. This drift is insignificant when compared with the total downrange distance of the 
trajectory. The cause of the drift is u nkn own and may be numerical. From a physical 
point of view, the aerodynamics of both the failed and nominal vehicles is symmetrical as 
is the model used for the trajectory. Closer examination of the control allocation and 
aerodynamic data may yield answers. One obvious difference between the two 
trajectories is that the nominal trajectory drifts slowly to the left while the vehicle with 
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stuck flaps experiences both left and right drift along the trajeetory. The angle of bank 
should be inspeeted to eorrelate the drift in erossrange. 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 20 Cross Range Profiles for Maximum Downrange Trajectories 

The altitude profdes have been diseussed in some detail previously, but an 
interesting event oeeurs near the end of eaeh trajeetory. Both trajeetories overshoot the 
endpoint altitude of 500 feet and appear to perform identieal maneuvers to meet the 500 
foot altitude requirement at the endpoint condition. This event was surprising not only 
because of the nearly identical maneuver, but also beeause the trajectory actually 
indicated a climb in the last moments of the trajectory. The encroachment below 500 
foot altitude was allowed by the problem formulation; the only limit plaeed upon the z 
state variable was that it had to be nonnegative. 
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Altitude, feet 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 21 Terminal Maneuver for Max Downrange Trajectories 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 22 Altitude Profiles for Max Downrange Trajectories 
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c. Velocity Profiles 

Figure 23 depicts the velocity profiles for the max downrange trajectories. 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 23 Velocity Profiles for Max Downrange Trajectories 


d. Flight Path and Azimuth Angle 

The flight path angles (FPA) shown in Figure 24 for these two trajectories 
differ as a consequence of the controls available to the individual vehicles. The failed 
vehicle exhibits a very steep descent but ends the trajectory with a terminal maneuver 
similar to the nominal trajectory. The similarity of the terminal maneuvers is a direct 
result of the endpoint conditions shared by the two trajectories. 
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Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



time, seconds 

Figure 24 Flight Path Angle Profiles for Max Downrange Trajectories 

There is no limit placed upon the azimuth angles for these trajectories and 
both vehicles wander somewhat in azimuth. The average azimuth of both trajectories in 
Figure 25 is positive and accounts for the left bias of the crossrange. Notice the nominal 
trajectory has the most stable azimuth between 100 and 700 seconds. The stable azimuth 
is small and positive and accounts for the steady left drift found in the nominal trajectory. 
The vehicle experiencing the stuck flap is far less stable and the azimuth varies by as 
much as 11 degrees. 
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Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 
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Azimuth Angle Profile for Max Downrange Trajectories 
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e. 


Vertical Speed Profile 


Nominal Trajectory 



Trajectory with Flaps Stuck at 26 degrees 



Figure 26 Vertical Speed Profiles for Nominal and Failed Trajectories (Max 

Downrange) 

f. Body Frame Normal Force 

The normal force experienced by the vehicle for both trajectories are 
depicted as functions of time and distance. As expected, the vehicle with the control 
surface failure has a more varied normal force profile with larger force excursions. Both 
vehicles have increases in normal force at the terminal maneuver as shown in Figure 27 . 
Additionally, the nominal vehicle is “pulling g’s” during the zoom climb maneuver early 
in the trajectory. 
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Symmetrical Stuck Flap Trajectory 80-nodes 



Figure 27 Normal Force (-1) vs. Time for Vehicle with Stuck Flaps 



20 40 60 80 100 120 140 
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Figure 28 Normal Force vs. Downrange Distance for Vehicle with Stuck 

56 








Flaps 


Nominal Trajectory 80-nodes 



Figure 29 Normal Force ( -1) vs. Time for Nominal Vehicle 


Nominal Trajectory 80-nodes 



Figure 30 Delta Normal Force( n^-l) vs. Downrange Distance for Nominal Vehicle 


57 










g. Controls 

While the problem formulation uses pseudo controls to induce inertia into 
the dynamical model, the physical controls of the vehicle are its angles of attack and 
bank. The angle of attack profiles clearly demonstrate the effects of the reconfigured 
aerodynamics and resulting angle of attack limits. The failed vehicle is limited to a much 
smaller range of angles of attack until the subsonic flight regime at which point it can use 
the same full range of angles of attack as the nominal trajectory. Both trajectories share 
the same terminal maneuver. The terminal maneuver is simply a climb that compensates 
for undershooting the 500 foot final altitude condition. This maneuver can be compared 
to the terminal maneuver shown in the altitude profile. 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 3 f Angle of Attack Profiles for Max Downrange Trajectories 

The bank angle profiles are as expected and compare well to the 
crossrange profiles previously discussed. The nominal vehicle is relatively steady all the 
way to the endpoint conditions where it gets a little chatter while the vehicle with the 
stuck flaps has bank angle excursions as large as 27 degrees. 
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Figure 32 Bank Angle Profiles for Max Downrange Trajeetories 


h. Pseudo Controls 

The pseudo eontrols depleted are not the true eontrols, but rather the 
eontrol rates. The profiles show that while the majority of the maneuvers request the 
eontrols at the maximum rate available (the pseudo eontrols are rate limited), not every 
maneuver requires full eontrol surface rates. The points where the pseudo controls are 
limited provide excellent insight into the optimality of the solution by examination of the 
necessary conditions. 
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Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 33 Angle of Attack Rate Profiles for Maximum Downrange Trajectories 



time, seconds 


Figure 34 Bank Angle Rate Profiles for Maximum Downrange Trajectories 
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L Trim Deficiency 

The restricted region in Mach-alpha space, due to the moment trim 
deficiency for the vehicle with stuck flaps, imposes severe trajectory performance 
penalties as previously shown. Figure 35 is meant to show the differences between the 
two trajectories in Mach-alpha space as well as to provide visual verification that the 
trajectory for the vehicle experiencing flap failure does not violate the imposed trim 
deficiency constraint of 0.003 and instead closely follows the contour in regions where 
the constraint is active. 


Trajectory Trim Deficiency 




Figure 35 Trim Deficiency Map for Max Downrange Trajectories 


j. Optimality Analysis 

The optimality of the two trajectories can be shown by investigating the 
Hamiltonian and necessary conditions. The Hamiltonian for these trajectories should be 
zero. Several problems inherent with the trajectory problem as posed can prevent DIDO 


61 






















from identifying a locally optimal solution and instead characterize the solution as “near 
optimal” where the cost may be within 99.999% of optimal. The two major problems 
which lead to “near optimal” solutions are 

1. The use of nonsmooth data. In this case, the aero data tables are not 
smooth and can cause some difficulty in the Jacobian matrices used in the DIDO code. 

2. Insufficient number of nodes used in the solution. This is a tradeoff 
between execution time and accuracy. If a 40-node solution yields propagation errors of 
acceptable magnitude along with a near-optimal solution, it may be undesirable to 
increase the execution time to find a locally optimal solution, if it exists. 

Both nominal and stuck flap trajectories used an 80-node solution which 
resulted in a near-optimal solution determination. The Hamiltonians are not exactly zero 
but are close. The stuck flap trajectory Hamiltonian is not as smooth as expected but 
investigation of the necessary conditions shows appropriate switching of the dual for 
constrained controls. 


80-node Nominal Trajectory, Max Down Range 



Figure 36 Hamiltonian for Max Downrange Nominal Trajectory 
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aoa rate, radians/sec Hamiltonian 


Symmetrical Stuck Flap Trajectory, 80-nodes 



Figure 37 Hamiltonian for Vehicle with (26 deg) Stuck Flaps 

Symmetrical Stuck Flap Trajectory, 80-nodes 



Figure 38 Constrained Pseudo Controls for (26 deg) Stuck Flap Trajectory 
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The pseudo controls for the vehicles experiencing the stuck flaps are 
generally constrained throughout the trajectory as shown in Figure 38 and Figure 39 The 
rapid switching between extreme pseudo control rates is a close approximation of bang- 
bang control. The necessary conditions for the pseudo controls are satisfied when the 
duals for both angle of attack and bank angle rates are positive for upper-bounded control 
rates and negative for lower-bounded control rates as shown in Figure 40 and Figure 41 


Symmetrical Stuck Flap Trajectory, 80-nodes 



Figure 39 Upper and Lower Bounded Bank Angle Rate for (26 degree) Stuck Flap 

Trajectory 
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controls dual 




Figure 40 Ra nk Angle Rate and Associated Dual for (26 deg) Stuck Flap Trajectory 


controls dual 



controls 



Figure 41 Angle of Attack Rate and Associated Dual for (26 degree) Stuck Flap 

Trajectory 
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k. Feasibility 

The feasibility of the trajectory is demonstrated by propagating the DIDO 
controls using an R-K algorithm and comparing the DIDO and propagated state vectors at 
the endpoint. There is an obvious tradeoff-a trajectory that seeks to place the vehicle in 
a landable attitude at the landing threshold of a runway needs sufficient accuracy at the 
endpoint conditions that does not place the vehicle in danger. The accuracy of an 80- 
node solution is generally sufficient. The errors in the corresponding components of the 
state vector are shown below. 

error = [436ft -267ft 79 ft 11.5 deg -10 deg -21 deg 0.11 deg/sec -0.63 deg/sec]’ 

A graphical depiction of the DIDO and propagated controls is shown in 

Figure 42 


DIDO and Propagated Controls for Stuck Flap Trajectory - 80 nodes 

1 6 I-1-^-1-^^^^^-r 



2 1_^^^^^^^^^_I 

0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

time, secs 

Figure 42 DIDO and Propagated Controls (DIDO controls are denoted by discrete 

points) 
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4, Minimum Downrange 

a. Mach-Altitude Profile 

The Mach-altitude profile for the minimum downrange trajectory 
produced some unexpected results. The nominal and stuck flap trajectories closely 
approximate each other until approximately 60,000 feet where the nominal vehicle 
conducts a zoom climb that approaches stall and then rapidly descends at a high subsonic 
speed. The point of the climb corresponds to altitude of minimum speed of sound in 
Figure 7 



Mach number 


Figure 43 Mach-altitude Profiles for Minimum Downrange Trajectory 
b. Trajectory in Cartesian Coordinates 

The profiles in Cartesian coordinates for the minimum downrange 
trajectory are shown in Figure 44 . The stuck flap trajectory shows better minimum 
downrange performance than the nominal trajectory. This is surprising because the 
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nominal vehicle has the full available range of control effector positions available while 
the failed vehicle only has a subset of the available control effector positions. In other 
words, the performance of the nominal vehicle should never be less than that of the 
vehicle experiencing a control surface failure. One explanation for the performance 
difference lies with the control allocation found in the inner loop of the algorithm. Recall 
that the control effector position vector for the requested angle of attack is optimized to 
minimize pitching moment using the aero dynamic data and vehicle flight condition 
(determined by Mach and angle of attack) as arguments. The resulting control effector 
position vectors are different with each aerodynamic configuration for the same requested 
angle of attack with accompanying differences in the lift and drag forces experienced by 
the vehicles. The consequence is nearly 10 miles in additional downrange distance and 
100 additional seconds of flight time. 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 44 Downrange Distance Profiles for Minimum Downrange Trajectories 

The crossrange profile of the vehicle with stuck flaps veers nearly a mile 
off centerline before attempting to return to centerline while the nominal trajectory 
remains on centerline throughout the trajectory. 
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Figure 45 Crossrange Profiles for Minimum Downrange Trajeetories 


The altitude profile for the nominal trajeetory shows the previously 
pointed out climb at 50,000 feet. Both trajectories exhibit similar terminal maneuvers as 
those meeting the maximum downrange criterion. The vehicles pass through 500 feet 
and initiate climbs to meet the endpoint conditions. 
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Altitude, feet 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 46 Altitude Profile for Minimum Downrange Trajeetories 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 47 Terminal Altitude Maneuver for Minimum Downrange Trajectories 
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c. Velocity Profile 

The vertical speed profiles are nearly identical until the transonic region at 
which point they diverge. This divergence is a consequence of the nominal trajectory’s 
zoom climb that occurs at 50,000 feet. The final 30 seconds’ velocity profile for both 
trajectories is very similar. 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 48 Velocity Profiles for Minimum Downrange Trajectories 
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Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 49 Terminal Velocity Profile for Minimum Downrange Trajectories 


d. Flight Path and Azimuth Angle 

The terminal flight path angle profile for both trajectories is very similar. 
The failed vehicle’s final azimuth performance is far better than the nominal vehicle with 
some attempt at ground track reversal in evidence. 
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Azimuth Angle, degrees Ciq’ Flight Path Angle, degrees 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 
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Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 
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Figure 51 Azimuth Angle Profile for Min Downrange Trajectories 
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e. 


Vertical Speed Profile 


Nominal Trajectory 



distance, miles 

Trajectory with Flaps Stuck at 26 degrees 



distance, miles 


Figure 52 Vertical Speed Profiles (vs. Downrange Distance) for Minimum 

Downrange Trajectories 
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Nominal Trajectory 
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Figure 53 Vertical Speed Profiles (vs. time) for Nominal and Failed Trajectories, 

Minimum Downrange 
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f. Body Frame Normal Force 

Normal Force for Stuck Flap Trajectory, 80 nodes 



Figure 54 Body Frame Normal Foree Stuek Flap Minimum Downrange 

Trajectory 
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Normal Force for Nominal Trajectory, 80 nodes 



Figure 55 Body Frame Normal Foree Nominal Minimum Downrange 

Trajectory 
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Controls 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 56 
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Figure 57 Bank Angle Profile, Minimum Downrange Trajeetories 
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Pseudo Controls 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 58 Angle of Attaek Rate Profile, Minimum Downrange Trajectories 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 59 Bank Angle Profile, Minimum Downrange Trajectories 
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Trim Deficiency 


Trajectory Trim Deficiency 




Figure 60 Trajectory Trim Deficiencies, Minimum Downrange Trajectories 


The trim deficiency path constraint for the failed vehicle is not active until 
passing Mach 2 at which point the vehicle follows the path constraint until subsonic. The 
initial Mach-alpha profiles for both vehicles are nearly identical until Mach 6 at which 
point they diverge. The vehicles’ Mach-altitude profiles begin to slightly diverge until 
Mach 4 at which point the nominal vehicle conducts a climb. While it appears to be 
some encroachment upon the trim deficiency path constraint for the failed trajectory 
between Mach 2 and Mach 1, all trim deficiency points corresponding to the trajectory’s 
nodal points are within limits. The possible encroachment illustrates the need to have 
sufficient nodes not only for endpoint accuracy but also to adequately characterize the 
trajectory as related to trim deficiency. The solution can avoid encroachment by placing 
knots at the critical points where there is a large change in trim deficiency for small 
changes in either Mach number or angle of attack. 
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j. Optimality Analysis 

The trajectories generated by DIDO were “near-optimal”. An examination 
of the Hamiltonians for both trajectories shows larger Hamiltonians than the expected 
zero, however, the average Hamiltonian for both trajectories is nearly zero. Both 
vehicles’ flight path angles reach the limiting value of ninety degrees. The duals for both 
FPAs meet the necessary conditions. 


Hamiltonian for Nominal Minimum Downrange Trajectory, 80 nodes 



Figure 61 Hamiltonian for Nominal Minimum Downrange Trajectory 
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flight path angleDual 



States 



Figure 62 Constrained Flight Path Angle and Associated Dual, Nominal Minimum 

Downrange Trajectory 


controls dual 



controls 



Figure 63 Constrained Angle of Attack Rate and Associated Dual, Nominal 

Minimum Downrange Trajectory 
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Scaled flight path angle Dual ^ Hamiltonian 


Hamiltonian for Stuck Flap Minimum Downrange Trajectory, 80 nodes 



64 Hamiltonian for Stuck Flap, Minimum Downrange Trajectory 


flight path angleDual 



States 



Figure 65 Constrained Flight Path Angle and Associated Dual, Stuck Flap Minimum 

Downrange Trajectory 
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controls dual 



controls 



Figure 66 Constrained Angle of Attack Rate and Associated Dual, Stuck Flap 

Minimum Downrange Trajectory 

k. Feasibility 

Both nominal and stuck flap trajectories propagate with acceptable error 
with 80-node solutions. Propagation errors are identical to those for the failed, maximum 
downrange trajectory. 
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DIDO and Propagated States for Stuck Flap Trajectory, 80 nodes 



time, TU 


Figure 67 DIDO and Propagated Controls for Stuek Flap, Minimum Downrange 
Trajectory (DIDO controls are denoted by discrete points) 

5, Minimum /Maximum Crossrange 

The minimum and maximum crossrange profiles show symmetry and as such 
only the maximum crossrange will be shown as representative trajectories. 

a. Mach-Altitude Profile 

The Mach-altitude profile is similar to those shown for the maximum 
downrange case. The vehicle trades airspeed for altitude as it seeks to maximize the cross 
range. 
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Mach - Altitude Profile for 80-Node Trajectory 



Figure 68 Downrange Distance Profiles for Maximum Crossrange Trajectories 
b. Trajectory in Cartesian Coordinates 

An interesting result of Figure 69 is that the nominal vehicle reaches the 
maximum downrange in the first fifth of its flight time. This corresponds to an azimuth 
of 90 degrees at the 300 second point. This failed vehicle initiates its turn at the same 
time as the nominal vehicle but is able to reach 90 degrees azimuth (with some trouble) 
much earlier because of the energy bleed caused by the additional drag of the stuck flaps. 


86 







Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 69 Downrange Distance Profile for Maximum Crossrange Trajectories 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 70 Crossrange Distance Profile for Maximum Crossrange Trajectories 
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Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 71 Altitude Profile for Maximum Crossrange Trajeetories 

c. Velocity Profile 

Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 72 Veloeity Profile for Maximum Crossrange Trajectories 
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d. Flight Path and Azimuth Angles 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 73 Flight Path Angle Profile for Maximum Crossrange Trajectories 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 74 Azimuth Angle Profile for Maximum Crossrange Trajectories 
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e. 


Vertical Speed Profiles 


Nominal Trajectory 



Trajectory with Flaps Stuck at 26 degrees 



Figure 75 Vertical Speed Profiles for Maximum Crossrange Trajectories 
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f. Normal Force Profiles 


Normal Force for Nominal Max Crossrange, 80-nodes 



Figure 76 Normal Force (“ 1) foi" Nominal Maximum Crossrange Trajectory 


Normal Force for Failed Max Crossrange, 80-nodes 



Figure 77 Normal Force( n^-l) for Stuck Flap, Maximum Crossrange Trajectory 
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Angle of Attack, degrees 


g- 


Controls 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 78 Angle of Attack Profile for Maximum Crossrange Trajectories 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 79 Bank Angle Profile for Maximum Crossrange Trajectories 
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h. Pseudo Controls 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 80 Angle of Attack Rate Profile for Maximum Crossrange Trajectories 


Comparison of Nominal and Stuck Flap Trajectories, 80-nodes 



Figure 81 Bank Angle Rate Profiles for Maximum Crossrange Trajectories 
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Mach 
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Trim Deficiency 


Trajectory Trim Deficiency 



Angle of Attack, Degrees 



Figure 82 Trim Deficiencies for Maximum Crossrange Trajectories 
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Hamiltonian 


7- 


Optimality Analysis 



Figure 83 Hamiltonian for Stuck Flap, Maximum Crossrange Trajectory 


Hamiltonian for 80 node Nominal Trajectory 



Figure 84 Hamiltonian for Nominal Maximum Crossrange Trajectory 
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Figure 85 Bank Angle Rate and Dual for Stuek Flap Max Crossrange Trajeetory 
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k. Feasibility 


DIDO and Propagated States for Norn Trajectory, 80 Nodes 
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Figure 86 DIDO and Propagated Controls for Stuck Flap Max Crossrange Trajectory 
(DIDO controls are denoted by discrete points) 


The error between the DIDO states and propagated states for the vehicle 
with stuck flaps is shown below. 


error = 


37/t 


-85/t -1ft -7 deg 9 deg -0.01^^^ -0.01 

/sec /sec / sec 


C. TRAJECTORY TO A SPECIFIED FIELD 

1, Endpoint Conditions and Cost Function 

The endpoint set for a trajectory to a specified field within the footprint contains 
the field Cartesian coordinates as well as the azimuth of the landing runway. Both 
nominal and failed vehicles were able to land at specified fields within their respective 
footprints however the choice of runway azimuths was limited to an approximate range 
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of +/- 100 degrees. For example, landing at a runway within the failed footprint with a 
runway heading of -90 degrees was aeeomplished relatively rapidly with a resultant 
“loeally optimal solution found”. The same field with a runway heading of -100 degrees 
gave a “near optimal” solution while a field with a runway heading of -120 degrees gave 
infeasible or suboptimal results for the same 80-node run. It should be noted that while 
the infeasibility was for an 80-node trajeetory, larger node trajeetories have been shown 
to provide near optimal or loeally optimal solutions when bootstrapped from lower-node 
infeasible trajeetories. There is the additional eomputation time penalty that eomes along 
with eomputing trajeetories with larger nodes. 

The eost funetion for these trajeetories to speeified field was a Lagrangian eost 
funetion that minimized the angles of attaek and bank. The rationale behind a eost 
funetion sueh as this is that depending on the eontrol system failure, exeessive movement 
of eontrol surfaee hydraulie aetuators may deplete the hydraulie working fluid of that 
partieular loop resulting in easeading failures. 

The full problem formulation is shown below: 

min„ T(x(t),M= j 

h) 

subjeet to x = f (x, u, r) 
hi - 

u,<u{t)<u^ 

Where 

= yf Zf 

= [110 miles -1 mile 500 feet -500±1000 ft/min Mach 0.15±2% -100° 0±10°]'' 
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All other initial conditions and path, state and control vectors and limits are 
identical to the previous formulations. 

2, Profiles for a Failure-Mode Trajectory to a Specified Field 
a. Mach-Altitude Profile 


Mach - Altitude Profile for 80-Node Trajectory 



Figure 87 Mach-altitude Profile for Stuck Flap Trajectory to Specified Field 
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Downrange Distance, miles 


b. Trajectory in Cartesian Coordinates 

80-node Trajectory to Specific Field 



Figure 88 Downrange Profile for Stuck Flap Trajectory to Specified Field (Flaps 

Stuck at 26 degrees) 
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80-node Trajectory to Specific Fieid 



Figure 89 Crossrange Profile for Stuck Flap Trajectory to Specified Field (Flaps 

Stuck at 26 degrees) 



time, seconds 


Figure 90 Altitude Profile for Stuck Flap Trajectory to Specified Field 
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3-Dimension Piot of Entry Trajectory 



Figure 91 Three Dimensional Trajeetory Depletion (not to seale) for Failed Vehiele 

(Flaps Stuek at 26 degrees) 
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Figure 92 Quiver Plot (not to scale) for Failed Vehicle Trajectory to Specified Field 

(Flaps Stuck at 26 degrees) 
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Ground Track Plot of Trajectory 



Figure 93 Ground Track Plot (not to scale) of Failed Vehicle Trajectory to Specified 

Field (Flaps Stuck at 26 degrees) 
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Velocity, mph 


C, 


Velocity Profile 


80-node Trajectory to Specific Field 



Figure 94 Velocity Profile for Failed Vehicle Trajectory to Specified Field 
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d. Flight Path and Azimuth Angle 


80-node Trajectory to Specific Fieid 



Figure 95 Flight Path Angle Profile for Failed Vehicle Trajectory to a Specified 


Field 

80-node Trajectory to Specific Fieid 



Figure 96 Azimuth Angle Profile for Failed Vehicle Trajectory to a Specified Field 
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e. 


Body Frame Normal Force 


stuck Body Flap, 80 node Trajectory 



Figure 97 Normal Force( n^-l) and Angle of Attack Profile for Failed Vehicle 

Trajectory to Specified Field 
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stuck Body Flap, 80 node Trajectory 



Figure 98 Normal Force( n^-l) and Angle of Bank Profile for Failed Vehicle 

Trajectory to Specified Field 
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f. Vertical Speed Profile 


Vertical Speed For 26 deg Stuck Flap Trajectory to Specific Field 



distance, miles 

Figure 99 Vertical Speed Profile for Failed Vehicle Trajectory to Specified Field 
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Vertical Speed For 26 deg Stuck Flap Trajectory to Specific Field 



Figure 100 Vertical Speed Profile (vs. time) for Failed Vehicle Trajectory to Specified 

Field 
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Controls 


g- 


80-node Trajectory to Specific Fieid 



Figure 101 Angle of Attaek Profile for Failed Vehiele Trajeetory to Speeified Field 


80-node Trajectory to Specific Fieid 



Figure 102 Angle of Bank Profile for Failed Vehicle Trajectory to Specified Field 
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Pseudo Controls 


80-node Trajectory to Specific Fieid 



Figure 103 Angle of Attack Rate Profile for Failed Vehicle Trajectory to Specified 

Field 
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Bank Angle Rate, degrees/second 


The bank angle rate profile is nearly bang-bang for this trajectory. 
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Figure 104 Bank Angle Rate for Failed Vehicle Trajectory to Specified Field 


113 
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Trim Deficiency 


Trajectory Trim Deficiency 




Figure 105 Trim Deficiency Profile for Failed Vehicle Trajectory to Specified Field 


The trim deficiency profile illustrates several problems associated with 
fitting curves to the trim deficiency data. The contour plot of the fitted trim deficiency 
curve is relatively well behaved except for the low Mach and angle of attack 
combinations, which is exactly where the trajectory wants to go. For this case, the 
process of fitting the trim deficiency data has artificially reduced the Mach-alpha 
envelope available to the trajectory with a corresponding decrease in the performance of 
the vehicle. There must be sufficient nodes or properly placed knots as previously 
discussed to cover regions in the trajectory where there is rapidly changing trim 
deficiencies. While the solutions points do not violate the trim deficiency path constraint 
as shown below, the trajectory between the points may violate the constraint. This 
violation can be easily removed with proper problem setup. 
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Columns 1 through 15 


0.0005 0.0004 0.0004 

0.0003 

0.0002 

0.0001 

0.0000 

0.0001 

0.0001 

0.0002 

0.0004 

0.0005 

0.0006 

0.0003 ■ 

0.0008 

Columns 16 through 30 
-0.0008 0.0002 0.0005 

0.0004 

0.0001 

-0.0002 

0.0001 

0.0014 

0.0025 

0.0030 

0.0030 

0.0030 

0.0024 

0.0030 

0.0030 

Columns 31 through 45 
0.0021 0.0022 0.0030 

0.0030 

0.0030 

0.0030 

-0.0020 

-0.0033 

-0.0013 

0.0011 

0.0030 

0.0030 

0.0011 

-0.0025 

-0.0081 

Columns 46 through 60 
-0.0145 -0.0182 -0.0201 

-0.0278 

-0.0412 

-0.0591 

-0.0806 

-0.1083 

-0.1281 

-0.1283 

-0.1106 

-0.0766 

-0.0514 

-0.0352 

: -0.0412 

Columns 61 through 75 
-0.0585 -0.0721 -0.0858 

-0.0912 

-0.0861 

-0.0600 

0.0002 

0.0001 

-0.0000 

-0.0001 

-0.0002 

-0.0001 

-0.0002 

-0.0001 

-0.0001 


Columns 76 through 80 
- 0.0000 0.0000 0.0001 0.0001 0.0001 

Table 2. Numerical Trim Deficiency Constraint Results for 80 nodal Points 


j. Optimality Analysis 


stuck Body Flap, 80 node Trajectory 



Figure 106 Flamiltonian for Failed Vehicle Trajectory to Specified Field (Flaps Stuck 

at 26 degrees) 
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Scaled flight path angle Dual 


Figure 107 


flight path angleDual 




FPA and Associated Dual for Failed Vehicle Trajectory to Specified Field 
(Flaps Stuck at 26 degrees) 
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Scaled bank angle rate 


controls dual 




Figure 108 Bank Angle Rate and Associated Dual for Failed Vehicle Trajectory to 
Specified Field (Flaps Stuck at 26 degrees) 
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Scaled aoa rate 


controls dual 




Figure 109 Angle of Attack Rate and Associated Dual for Failed Vehicle Trajectory 
to Specified Field (Flaps Stuck at 26 degrees) 
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k. Feasibility Analysis 


DIDO and Propagated States for Stuck Flap Trajectory, 80 Nodes 



time,TU 


Figure 110 Scaled DIDO and Propagated States, Failed Vehicle Trajectory to 
Specified Field (DIDO controls are denoted by discrete points) 


The error between the DIDO states and the R-K propagated states using 
the DIDO controls is shown in Figure XX. The altitude error is still high and either the 
trajectory should be generated using more nodes or the endpoint conditions relaxed 

Error = [108ft 1247ft -460ft -01.7ft/sec -O.Sdeg 4.5deg Odeg/sec -O.ldeg/sec] 
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D, DIDO PERFORMANCE 

While there are several performance criteria for the optimal trajectory solutions, 
several criteria dominate; 

i. Accuracy of the solution. In other words, how closely do the states at the 
endpoint set generated by using the propagated solution controls match the 
internally generated states used by the optimization algorithm. 

ii. Optimality. How optimal is the solution? Examination of the necessary 
conditions can help answer this question but can be inadequate when 
attempting to qualify a “near optimal” solution. 

iii. Time to find a solution. Can the solution be found in “real time” and if 
not, how close to real time can a solution be found? 

For the psuedospectral method used in this thesis, all three criteria are a function 
of the number of nodes used in the solution. Given a feasible solution, the error between 
propagated and DIDO states can be reduced by increasing the number of nodes. At the 
same time, increasing the number of nodes in the solution comes at significant expense in 
computation time. If a locally optimal solution is declared, the end states will not 
generally change; however, if a “near optimal” solution is declared, then the end states 
become a function of the number of nodes in the solution. Recall that the generation of a 
“near optimal” solution is a consequence of many things including insufficient number of 
nodes and the use of nonsmooth data. 

The following analysis provides some insight into the DIDO solution performance 
for a nominal trajectory with the same initial and endpoint conditions as shown in the 
results. The graphs represent data from a nominal maximum downrange trajectory. The 
relationship between propagated errors and number of nodes in the solution is shown in 
Figure III and Figure 112 Notice that the error for low node solutions is unacceptable 
while the error for solutions greater than 60 nodes results is acceptable. 
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error, miles 


Error Between DIDO and Propagated States 



Figure 111 Downrange and Crossrange Error, Nominal Trajeetory 


Error Between DIDO and Propagated States 



Figure 112 Altitude Error vs. Nodes 
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DIDO execution times for P4 2.0 Ghz Processor 



Figure 113 Time to Reach a Solution, Nominal Trajectory for Pentium 4, 2.0 Ghz 

processor using Windows 2000. 

While larger nodes give better propagation performance, the cost is significant. 
The execution time for an 80-node solution is nearly double that of a 60-node solution. 
These times were for unoptimized code running on a networked Pentium 4 processor 
running Matlab 6.0 in a full Windows 2000 environment. Considerable performance 
increases resulting from speed-optimized code and a faster processor using a dedicated 
(non-multitasking) computing environment can be expected. Execution times of 70 
seconds for 20- node nominal trajectories have been routinely observed on a P4 3.06 Ghz 
laptop with minimum cpu overhead. 

Figure 114 and Figure 115 illustrate the convergence of the “near optimaf’ 
solution. Notice that by 60-nodes the final DIDO cross and downrange state variables are 
nearly identical while altitude convergence is reached immediately. The intent of these 
graphs is to show the pitfall of choosing a small number nodes to generate faster 
solutions. The difference between the maximum downrange of a 20 node solution and a 
60 node solution is nearly ten miles. 
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Convergence of Near Optimal Solutions 



Figure 114 Convergence of Down and Cross Ranges for Near Optimal Solutions 
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V. FUTURE WORK 


A, CONVERSION TO SI UNITS 

This thesis was conducted using the English unit system. The units for velocity, 
heating rate and pressure do not lend themselves to convenient balancing or scaling. 
Ross [ref 21] has shown that the use of English units in the problem formulation can 
result in numerical difficulties leading to infeasible solutions. The metric system is a 
more intuitive system with all units a subset of the primary distance, mass and time units 
of meters, kilograms and seconds. Problems formulated using the inks system lead to 
more easily balanced and scaled units with better corresponding numerical results. 


B, DATA SMOOTHING 

The DIDO software package performance is best with smoothed data that avoids 
singularities in the gradients of the state variables. The presence of non smooth data 
results in extremely long calculation times, “near optimal” vice “locally optimal” 
solutions along with the presence of “sweet spots” where a small range of nodes results in 
feasible solutions. One way to artificially smooth the data is to increase the number of 
nodes in the solution; however, the computation penalty can be substantial. The best 
approach may be to use a more finely meshed lookup Cd and Cl table in the dynamics 
file. These tables are generated offline and have minor impact on the lookup routine; the 
extra lookup time is easily offset by the quicker convergence time. 


C. TRIM DEFICIENCY MAPPING 

The trim deficiency function used in this thesis was a polynomial fitted to the 
tabular trim deficiency data. This was done to minimize execution time of the solution. 
The curve fit to the data was not perfect and a better approach may be to use a finely 
meshed lookup table as discussed in the previous section. This would help automate the 
trim deficiency process with some penalty in execution time. 


125 



D, CLOSED LOOP IMPLEMENTATION 

While 80 node solutions can take up to 50 minutes to converge on a Pentium-4 
2.0 ghz processor, an optimized code using 40 nodes can achieve convergence in as little 
as 1 minute on the same machine. While 40-node solutions do not generally provide the 
necessary accuracy for the endpoints, it may be possible to reduce outer loop execution 
time using lower order solutions and still meet the accuracy requirements by using 
feedback. This should be more closely investigated to determine the feasibility of this 
approach. 


E. OPTIMIZING THE DIDO CODE FOR SPEED 

The algorithms used in this thesis were constructed with minimal emphasis on 
computation speed. There are many options available to increase the speed of the 
algorithm. The first and most significant speed increase can be seen by removing the 
code from the Windows (multi-tasking) environment and compiling for use on a 
dedicated real time processor. Additional improvements in the MATLAB code can 
provide further speed increases; MATLAB 7 shows some speed increase over MATLAB 
6.5. As CPU speed continues to increase and the operating code is optimized, execution 
times should continue to decrease. Additional future improvements to the DIDO 
software package have the potential of achieving significant decreases in solution times. 


F, FULL 6-DOF MODEL WITH INNER AND OUTER LOOP COUPLING 

As detailed in the footprint discussion, the decoupled inner and outer loops along 
with a reduced order model do a satisfactory job of computing the vehicle performance. 
The true performance of the vehicle; however, is not realized because the two loops are 
decoupled as illustrated in the differences in the nominal and failed vehicle trajectory 
footprints. While the use of a 6-DOF model with inner and outer loop may appear to be 
the more difficult problem to solve, the complex interplay between the dynamical 
equations and state and path constraints may reduce the problem to one more easily 
solved. The potential benefits of having a fully integrated inner and outer loop are great 
enough to warrant further investigation into the feasibility of this approach. 
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